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Abstract 

A modified formulation of the Schrodinger functional (SF) is proposed. In the con- 
tinuum it is related to the standard SF by a non-singlet chiral field rotation and 
therefore referred to as the chirally rotated SF (%SF). On the lattice with Wilson 
fermions the relation is not exact, suggesting some interesting tests of universality. 
The main advantage of the xSF consists in its compatibility with the mechanism 
of automatic 0(a) improvement. In this paper the basic set-up is introduced and 
discussed. Chirally rotated SF boundary conditions are implemented on the lattice 
using an orbifold construction. The lattice symmetries imply a list of counterterms, 
which determine how the action and the basic fermionic two-point functions are 
renormalised and 0(a) improved. As with the standard SF, a logarithmically diver- 
gent boundary counterterm leads to a multiplicative renormalisation of the fermionic 
boundary fields. In addition, a finite dimension 3 boundary counterterm must be 
tuned in order to preserve the chirally rotated boundary conditions in the interacting 
theory. Once this is achieved, 0(a) effects originating from the bulk action or from 
insertions of composite operators in the bulk can be avoided by the mechanism of 
automatic 0(a) improvement. The remaining 0(a) effects arise from the boundaries 
and can be cancelled by tuning a couple of 0(a) boundary counterterms. The gen- 
eral results are illustrated in the free theory where the Sheikholeslami-Wohlert term 
is shown to affect correlation functions only at 0(a 2 ), irrespective of its coefficient. 
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1 Introduction 



The Schrodinger functional (SF) [IH3] has become a general tool to address non- 
perturbative renormalization problems in lattice QCD [4^*1 and is now even used 
in candidate theories which might describe a strongly interacting electroweak sec- 
tor |13H17j . Renormalization schemes based on the SF are gauge invariant, quark 
mass independent through renormalization in the chiral limit and suitable for eval- 
uation by both Monte Carlo and perturbative methods. Moreover, as the finite 
space-time volume is used to set the renormalization scale, recursive finite size tech- 
niques can be applied to bridge large scale differences, thereby avoiding to resolve 
widely different scales on a single lattice. As systematic errors in large scale lattice 
simulations are becoming smaller, the need for a better controlled continuum limit 
will likely stir more interest in SF schemes, and it is therefore necessary to address 
its remaining shortcomings. 

A first problem consists in the difficulty to implement the SF boundary con- 
ditions for fermion regularisations other than Wilson fermions. This has been 
largely solved, and formulations of the SF with staggered and Ginsparg- Wilson-type 
fermions have been given in |18H20j and |21H25| . respectively. 

The main problem with SF schemes is the presence of lattice artefacts which 
are linear in the lattice spacing a. Some of these 0(a) effects are caused by the mere 
presence of the Euclidean time boundaries, together with local Dirichlet conditions 
for the fields. Such 0(a) boundary effects will be present with any regularisation, 
since they do not arise from the breaking of a continuum symmetry. In practical 
applications, these effects are cancelled by adding a couple of boundary counterterms 
to the action, with coefficients determined perturbatively up to two-loop order [26\ 
127] . While this is often sufficient in practice, a non-perturbative determination of 
such boundary counterterms would be desirable and is indeed conceivable. 

With Wilson quarks there is a second category of 0(a) effects which are can- 
celled by the usual 0(a) counterterms to the Wilson quark action and the composite 
fields which appear in the correlation functions. At first sight it may be surprising 
that these standard 0(a) counterterms are required at all, given that massless Wil- 
son quarks in a finite volume enjoy the property of automatic 0(a) improvement [28]. 
However, the argument for automatic 0(a) improvement relies on a discrete chiral 
symmetry which is expected to be recovered in the continuum limit. The argu- 
ment fails as the standard SF boundary conditions break chiral symmetry, so that 
observables cannot be classified as either even or odd under this symmetry |29|30| . 

In this paper a modified definition of the Schrodinger functional for Wilson- 
type quarks is introduced, which is suitable for QCD with an even number of quark 

1 see [5] for early references and [6til2| for a selection of more recent applications and further 
references. 
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flavours, and which is compatible with automatic 0(a) improvement. In the contin- 
uum limit this modified SF is related to the standard SF by a chiral field rotation 
and is therefore referred to as the chirally rotated SF (%SF). A first account of this 
work has been given a while ago in refs. p9 30J. The paper is organised as follows: in 
Section 2 the argument of automatic 0(a) improvement is reviewed and the reason 
why it fails in the presence of SF boundary conditions. Possible modifications of 
the SF boundary conditions which may restore this argument are discussed next. 
Concentrating on the solution provided by the chirally rotated SF, its lattice imple- 
mentation through an orbifold construction is explained in Section 4. The lattice 
symmetries imply the counterterm structure relevant for renormalisation and 0(a) 
improvement (Section 5). In Section 6, the renormalisation constants and boundary 
0(a) improvement coefficients are determined to tree-level of perturbation theory, 
and automatic bulk 0(a) improvement is demonstrated to this order. The paper 
ends with some concluding remarks and an outlook (Section 7). Finally, Appendix A 
demonstrates the absence of zero modes in the continuum Dirac operator and Ap- 
pendix B collects some useful formulae for the free fermion propagator both on and 
off the lattice. 



2 Automatic O(a) improvement and SF boundary conditions 

In order to understand how the SF boundary conditions interfere with 0(a) improve- 
ment, it is useful to recall the argument why massless Wilson fermions are expected 
to be automatically 0(a) improved in a finite volume without boundaries. 

2.1 Symanzik's effective theory and automatic O(a) improvement 

The Wilson fermion action is given by 

S f = a^^xXAv + moM*), (2-1) 

X 

3 . 3 

Av = EH^ + ^hc-iW + ^o °^F^ x )- ( 2 - 2 ) 

/j,=0 fi,i>=0 

Here, niQ is a bare mass parameter, the covariant lattice derivatives in the Wilson- 
Dirac operator are defined as usual, and the last term is the Sheikholeslami-Wohlert 
term with F^ u denoting the clover leaf definition of the field strength tensor (see 
ref. [3T] for conventions and unexplained notation). Standard Wilson quarks are 
obtained for c sw = 0. We will focus on the massless limit, which is obtained by 
tuning the bare mass parameter to its critical value, niQ = m CT . 
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Cutoff effects can be studied using Symanzik's effective continuum theory 
[31] . The effective action, 

S cS = S + aSx + a 2 S 2 + 0(a 3 ), (2.3) 
starts out with the massless QCD continuum action, 

So = "27 J d * x te{FA x ) F A x )} + J ^ x $(xm(x). (2.4) 

The terms which may appear in Si must be integrals over local composite fields of 
mass dimension 5. In the chiral limit there are only two such terms, 



Si = ici J d x ip(x)a fJiU F^(x)ip(x) + c 2 J d x ip(x)D^ l D fJi ^(x), (2.5) 

which can be reduced to a single term by restricting attention to on-shell quantities. 
The effective action renders the cutoff dependence explicit, up to logarithmic terms 
(cf. ref. |33j). The leading cutoff effects in a renormalised connected correlation 
function of a (multi-local) observable O are then described by the corresponding 
continuum correlation functions with insertions of either Si or the 0(a) counterterm 
50 for the observable O, 

(0) c = (0}™ nt -a(SiO)™ nt + a(50) c c ont + 0(a 2 ). (2.6) 

Here, (-)c° nt refers to Euclidean expectation values taken with respect to the contin- 
uum action So, only taking into account connected diagrams. Obviously, the contin- 
uum action, So, is chirally symmetric. In particular, the discrete 75-transformation 

rp -4 7 5 V>, $ -)• -$75, (2.7) 

leaves So invariant, while it changes the sign of Si. Since the 75-transformation 
squares to the identity, one may decompose any composite field into parts which are 
either even or odd under this transformation. Note that, despite the loss of chiral 
invariance on the lattice, such a decomposition can be made unambiguously for any 
lattice field. An important observation then is that the 0(a) counterterm, 50 even 
for any 75-even field O ev en must be 75-odd and vice versa [28]. Provided that the 
functional measure is invariant under the change of variables to the 75-transformed 
fields, one is then led to conclude that 

(£lO cv cn)c° nt = — (SlOevcn)c° nt = ^' (^) 

(SO cven ) c c ont = -(5O even )l oat = 0, (2.9) 
i.e. Eq. (12. 6p for 75-even observables reads 

(Oevcn)c = (Oeven)™' 11 + 0(« 2 ). (2.10) 
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This does not mean that 0(a) effects are absent, rather these are present in cor- 
relation functions of 75-odd observables 0( jdi for which the leading term in (|2.6|) 
vanishes, so that 

(Oodd)c = -a(SiOod d r nt + a{50 odd ) c c ont + 0(a 3 ). (2.11) 

Note that the corrections here are of order a 3 , rather than a 2 . In fact, the mechanism 
of automatic 0(a) improvement more generally implies that lattice effects in 75- 
even (75-odd) correlation functions only come with even (odd) powers of the lattice 
spacing a |28|34| . 

It should be emphasised that the inherent 0(a) ambiguity in the definition of 
the chiral limit with Wilson quarks does not alter this picture. To understand this, 
suppose that mo has been tuned critical e.g. by requiring the divergence of the axial 
current to vanish for some particular matrix element. By making a different choice, a 
non-zero quark mass m' = aA 2 is generated relative to the previous definition, with 
A denoting some mass scale. Its effect can be accounted for by an insertion of a 
fermionic mass term. Since this term is 75-odd, its insertion into an even correlation 
function vanishes, 

J d 4 x (O even 4;(x)4,(x)) c ; nt = 0. (2.12) 

In order to obtain a finite contribution to 75-even observables at least a double 
insertion of this term is required, making it an a 2 -effect, as anticipated. 

Note that the above arguments rely on the assumption that the functional 
measure is invariant under the change of variables to the 75 -transformed fields. 
To check this assumption it is useful to think of the effective Symanzik theory as 
being formulated on a very fine lattice with a spacing b very much smaller than 
a such that the theory is essentially continuum-like. Nevertheless, the functional 
integral in a finite volume is still finite-dimensional and mathematically well-defined. 
Furthermore, with Ginsparg- Wilson quarks, the fermionic part of the continuum like 
action 

S 0J = b 4 ^(x)D N ^(x), (2.13) 

X 

has an exact chiral symmetry [35] . Here, Dn is the Neuberger operator [36], which 
satisfies the Ginsparg- Wilson relation [37] . 

j 5 D N + D N75 = bD Nl5 D N . (2.14) 

Combined with the 75-conjugation property, D* N = 75-^^75, this gives rise to the 
definition of 75 , 

75 = 75^, V = 1 — bD N , (2.15) 



4 



which is hermitian and squares to the identity. However, in contrast to the ordinary 
75, it depends implicitly on the gauge field through the unitary matrix V. Replacing 
the 75-transformation (I2.7P with its version on the fine lattice, 

rp ->• 7 5 V>, 4> -)• -$75, (2.16) 

one readily verifies the invariance of the action (|2.13p . which follows from the 
Ginsparg- Wilson relation ()2.14p . Changing variables in the functional integral to 
the 75-transformed fields, one finds that the Jacobian J is given by 

J = det( 7 57 5 ) = det(V) = det(7 5 y7 5 ) = det(0) = ±1. (2.17) 

For Nf quark flavours, each 75-transformed flavour thus contributes a factor J, so 
that the total Jacobian becomes J f . Hence, for even Nf, the 75-transformation 
(|2.16p is an exact symmetry of the path integral measur^l. In conclusion, one 
expects that lattice QCD with an even number of Wilson type quarks is indeed 
automatically 0(a) improved, at least in a finite volume with periodic boundary 
conditions. 

However, the situation changes in the presence of SF boundary conditions, 
which, for the fermionic fields, read [2] , 

P + V(x)U =0 = 0, P^(x)\ X0=T = 0, 

$(x)P-\ Xo=Q = 0, ^{x)P + \ X0=T = 0, (2.18) 

with P± = ^(1 i 70) ■ The 75-transformation is then no longer a symmetry of the 
effective continuum theory, as the 75-transformed fields satisfy the SF boundary 
conditions with the complementary projectors, e.g. 

^' = 75^ P-^\x =o = 75-P+V>L=o = 0. (2.19) 

This may be taken as a property of the fermion measure, since the function space 
integrated over is not the same before and after the change of variables. Hence, the 
transformed correlation functions cannot be proportional to the original ones, and 
the basic argument needed for automatic 0(a) improvement cannot be made. 

2.2 Rendering automatic O(a) improvement compatible with the SF 

The question then arises whether there are ways to save automatic 0(a) improve- 
ment in the presence of SF boundary conditions. One needs to either find an al- 
ternative to the 75-transformation or modify the SF boundary conditions such that 

2 For a single quark flavour, a gauge field dependent sign in the Jacobian remains, thereby spoiling 
the argument for automatic 0(a) improvement. While some additional argument would be needed 
for Nf = 1, the situation is slightly better for odd flavour numbers Nt > 1, as the argument can 
then be applied to any pair of fermion flavours. 
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the transformation remains a symmetry of the massless continuum action, while the 
(homogeneous) SF boundary conditions are left invariant. There are various ways 
to achieve this: 

• Let ip be a flavour doublet, and augment the 75-transformation by a flavour 
permutation, i.e. 

if) -»■ t 1 ^, i> -»■ (2.20) 

where r 1 is a Pauli matrix. This flavour permutation does not affect the bulk 
action and 0(a) improvement for 75T 1 -even observables in the bulk can be 
shown as before. Then, giving the boundary conditions a flavour structure 
through the replacement 

P± = i(l±7o) — >• P ± = i(l±7or 3 ), (2.21) 

one easily sees that, due to 

[T8r\P±]=0, (2.22) 
the boundary conditions are preserved by the 75 r 1 -transformation. 

• Another possibility consists in combining the 75T 1 -transformation (|2.20p with 
Dirichlet boundary conditions using the projectors 

g ± = i(l±i 7 07 5 r 3 ) => fo i T 1 ,Q±] = Q. (2.23) 

It is this option which will be investigated in much more detail below. 

While it is relatively easy to find alternative SF boundary conditions and/or alter- 
natives to the 75-transformation, the difficult part is to show that this leads to a 
sensible alternative definition of the Schrodinger functional in the continuum limit. 
Besides renormalizability one would like to maintain the properties of the standard 
SF such as the absence of fermionic zero modes, which is crucial for practical ap- 
plications. Another difficulty lies in the actual implementation of these boundary 
conditions on the lattice. All these questions will be addressed in the following for 
the proposed solution involving the projectors in Eq. (|2.23p . 

3 Continuum considerations 

The boundary conditions involving the projectors Q± can be thought of as a chi- 
rally rotated version of the standard SF boundary conditions. As the non-singlet 
chiral rotation is part of the continuum symmetries of massless QCD, the alterna- 
tive formulation of the SF to be introduced here will in fact be closely related to the 
standard SF. It is useful to have a closer look at this relationship before passing to 
the lattice formulation. 
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3.1 Chirally rotated SF boundary conditions 

Let ip' and t/v be isospin doublets of quark and anti-quark fields satisfying the ho- 
mogeneous SF boundary conditions, Eq. (]2.18p . Performing a non-singlet chiral field 
rotation, 

ip' = R(a)tp, ip' = $R(a), R(a) = exp(iaj 5 r 3 /2), (3.1) 

the rotated fields satisfy the chirally rotated boundary conditions, 

P + (a)ip(x)\ Xo=0 = 0, P-(a)ip(x)\ Xo = T = 0, 

VKx) 7o P-(a)Uo=o = 0, $(x) l0 P + (a)\ X0=T = 0, (3.2) 

with the projectors 

P±(at) = \ [1 ± 70 exp(ia 75 r 3 )] . (3.3) 

For a = the standard projectors P± = -P±(0) are recovered, while for a = n/2 one 
finds, 

P±(ir/2) = Q± = ^(1 ± i7075T 3 ). (3.4) 

Since 7oQ± = Q^Jo, the boundary conditions (|3.2p in this special case read 

Q+ip(x) \ xo= o = 0, Q_^(x)| X0=T = 0, 

^(a;)Q+| xo =o = 0, v5(x)Q_| X0=T = 0. (3.5) 

We thus see that boundary conditions involving the projectors Q± appear naturally 
when chirally rotating quark and anti-quark fields which satisfy the standard SF 
boundary conditions. 

3.2 Properties of the continuum Dirac operator 

In the continuum the rotation angle a parameterizes a family of linear spaces of 
functions which satisfy homogeneous Dirichlet boundary conditions in Euclidean 
time. In the spatial directions we assume L-periodicity and we allow for a constant 
U(l) gauge background field i9^/L. Formally this field may be included in the Dirac 
operator by promoting the SU(iV) colour gauge field to a U(iV) gauge field as follows: 

A,(x) = Y, A a ,(x)T a , A°(s)T° = ifl. (3.6) 

a=0 

Following |38|31j . we set 9\ = 9 2 = 0% = and #0 = 0- With respect to the natural 
inner product, 

(<P,x) = J dx J d 3 xv?(x)t x (x), (3.7) 
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the Dirac operator for a quark doublet and generic angle a is 75-hermitian up to a 
flavour exchange, i.e. 

(75T 1 [Ip + m + ^ q75 T 3 ] ) f = 75T 1 [Ip + m + «/u q 7 5 r 3 ] . (3.8) 

For later reference both a standard and a twisted mass term have been included. 
Hence, a well-defined eigenvalue problem is obtained for any valu^f] of a. Denoting 
the eigenfunctions in the standard SF by v n , the eigenvalue equation reads 

75T 1 [Jp + m + 2/iq 7 5T 3 ] v n = X n v n . (3.9) 

Performing the chiral rotation as in Eq. (|3.ip . 

V n = R(a)w n 

and using 

R{a) [ip + m + ifi q ~/ 5 T 3 ] R(a) = Ip + fh + ifl^r 3 , (3.11) 

one finds, 

75T 1 [Ip + fh + i/I q7 5T 3 ] W n = \ n W n . (3.12) 

Hence, the eigenfunctions are related by a chiral rotation, and the eigenvalues are 
exactly the same as in the standard SF, provided the mass parameters either vanish 
or are transformed covariantly, according to 

fh = m cos a — fi g sin a, /I q = msina + fi q cosa. (3.13) 

Moreover, the absence of zero modes in the massless theory can be established 
following Liischer [24J, and some details have been deferred to Appendix A. Given the 
absence of zero modes, the free continuum quark propagator can then be calculated 
using standard techniques. Some results are collected in Appendix B. 

3.3 Relating correlation functions at different values of a 

The chiral rotation (|3.ip is not anomalous, so that one may perform a corresponding 
change of variables in the functional integral. One obtains relations between cor- 
relation functions defined with the standard and the chirally rotated SF boundary 
conditions. This is similar to and generalises the relations between correlation func- 
tions in twisted mass and standard QCD (39 30J. In particular, one expects that 
these formal continuum relations will hold up to cutoff effects, for appropriately 
renormalised correlation functions on the lattice. 

3 As an aside note that for a = ty/2 the massless operator ilp is hermitian for a single fermion 
with boundary conditions corresponding to either the up or the down flavour in Eq. (|3.5|) . 
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It is instructive to include both a standard and a twisted mass term in the 
fermionic action at a = 0, 

S f [m,fi q ]= [ dxo/ d 3 xft(x)(ip + m + ifi^ 5 T 3 )ij'(x), (3.14) 
io Jo 

even though the focus will later be on massless QCD. Performing the change of 
variables, Eq. (|3.1|) . in the Euclidean functional integral one obtains the generic 
identities, 

<°[^]W q ,P ± ) = W^W) (W± ( 8) ) . (3-15) 

where the integration variables are assumed to be ip and ijj on both sides of the 
equation. In addition to the boundary conditions for tp at xq = 0, T, the subscript 
to the correlation functions indicates the mass parameters in the action and the mass 
parameters are related as in Eq. f)3. 13|) . The relation (|3,15p establishes a mapping 
of correlation functions defined at different values of a. An equivalent form of this 
mapping is given by 

W'fl We ± ) = < [ B (- a )^'^(- a )]>(^ qi p ± ( f - a )) • ( 3 - 16 ) 

Setting the twisted mass term to zero in the chirally rotated SF, and choosing 
a = vr/2, this simplifies to 

<°[^]>foo,0 ± ) = <O[«(-T/2)V > ^(-T/2)]>(0,-m ) p ± ), (3-17) 

which shows that the standard mass parameter in the chirally rotated SF is equiv- 
alent to a (negative) twisted mass parameter in the standard SF. 

The simplest fermionic correlation function is the quark propagator. Using the 
explicit expressions of Appendix B, one may easily verify the identity (|3.ip in the 
free theory, viz. 

w*>K^ = R{a) <^)^)>(SU(«)) ^ 

Gauge invariant correlation functions corresponding to the ones defined in the stan- 
dard SF |31|40j can be defined as well. For this we need the fermionic boundary 
fields, C)C and CjC'j which are naturally included in the mapping f)3. 15|) by iden- 
tifying them with the non-Dirichlet components near the time boundaries. More 
precisely we set 

C(x) = P_(a)V(0,x), C'(x) = P + (a)^(T,x), 

C(x) = ^(0,x)P + (-a), C'(x) = ^(T,x)P_(-a), (3.19) 

thus leaving the a-dependence implicit. This should not lead to any confusion, as 
the fermionic boundary fields mostly appear as insertions into correlation functions 
with the subscripts indicating the projectors used for the boundary conditions. 
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3.4 Symmetries and conventions 

The standard SF boundary conditions break all chiral symmetries, as is best illus- 
trated by the Ward identities [21]. With the global axial and vector non-singlet 
variations, 

= 75 y V, ^ = yV>, 

= ^75y , = , (3-20) 

the Ward identities take the form, 

<*xO>(p±) = ( S * SO )(p±) > X = A ' V ' ( 3 ' 21 ) 

where O denotes an arbitrary composite field, and the subscript to the expectation 
values indicates the boundary conditions for ifi at xq = 0, T. The variation of the 
action is given by 

SxS = - J ^xd^x), (3.22) 

with the Noether currents 

Al{x) = ^(x) 7 ^75y ip(x), V£{x) = $(x) 7 ^(x). (3.23) 

The boundary conditions in the spatial directions being periodic, the integral over 
the total divergence yields only surface terms in the time direction. The r.h.s. of 
Eq. (|3~2D becomes 

(S^SO) (P±) = I d 3 x([X a (0,x) - X a (T,x)] 0) {P±) . (3.24) 

Using the standard SF boundary conditions, Eqs. (|2.18p and the identity 7o-P± = 
±P±, one obtains 

A a (0,x) = ^(0,x)P +75 yP-VK0,x) = C(x) 75 yC(x), (3.25) 
A a (T,x) = V5(r,x)P_ 75 yP + ^(0,x) = C , (x)75yC , (x), (3.26) 
whereas Vo a (0,x) and Vq(T, x) vanish. One thus arrives at the Ward identities, 

(KO) { p ±) = i|d 3 z<[C(z) 75 r a C(z) + C / ( Z )75r a C / (z)]0) (P±) , (3.27) 
<<5 v O> (P±) = 0, (3.28) 
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which show that the standard SU(2) vector symmetry is conserved, while chiral sym- 
metry is broken by terms which are localized at the time boundaries, with coefficients 



Proceeding analogously we may derive the Ward identities in the presence of 
chirally rotated SF boundary conditions, 



Using the projector structure implicit in the fermionic boundary fields, Eq. (|3.19j) . 
the r.h.s. of the vector Ward identity (|3.30p with index a = 3 is seen to vanish. 
Moreover, there are two further linear combinations of Eqs. (|3.29l3.30p for which 
the r.h.s. vanishes, so that three generators correspond to conserved symmetries. 
This is not surprising, as the chiral rotation merely transforms the vector symmetry 
of the standard SF to an unusual form, a situation which is familiar from twisted 
mass QCD |39|30| . To illustrate this further, consider the case of most interest, 
a = 7r/2, where the Ward identities take the form 



Apparently a mixture of vector and axial symmetries are broken by boundary terms. 
However, in massless QCD this is really a convention dependent statement. A 
preferred or "physical" basis is only available in the presence of a mass term: those 
symmetry transformations which leave the mass term invariant are identified as 
vector symmetries whereas the remaining chiral flavour transformations are classified 
as axial symmetries. In massless QCD with SF type boundary conditions we will 
choose a similar convention, and define the "physical" basis as the one with standard 
SF boundary conditions. This implies that the broken symmetries are in all cases 
the axial symmetries, whereas the vector or isospin symmetries are conserved. The 
situation is then completely analogous to twisted mass QCD, in particular, the 
"dictionary" translating composite fields between both bases can be taken over from 
twisted mass QCD. 

4 xSF boundary conditions on the lattice 

Implementing boundary conditions on the lattice is not straightforward, since bound- 
ary conditions cannot really be imposed. Rather, they arise dynamically from the 



of 0(1). 
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structure of the lattice action near the boundaries. An attractive, though not uni- 
versally applicable technique to obtain a lattice theory with the desired boundary 
conditions is based on an orbifold construction. In the context of lattice QCD orb- 
ifold techniques have previously been applied by Taniguchi [21] in order to implement 
SF boundary conditions for Ginsparg- Wilson quarks. The role of the orbifold con- 
struction is to define the action of the lattice operator near the time boundaries, 
such that the desired boundary conditions are obtained for the free theory. In a 
second step one needs to analyse the symmetries and list the allowed boundary 
counterterms which will be generated by the interactions. From the list of possible 
counterterms one may conclude whether the desired boundary conditions will be 
obeyed in the renormalised theory, or whether this requires the fine tuning of some 
counterterm coefficients (cf. Sect. 5). 

4.1 Orbifold construction 

The starting point is a single Wilson quark flavour x with lattice action 

S f [x,X,U] = a 4 £ ^x{x){Dw + mo) X (x). (4.1) 

-T<x <T x 

The fermionic fields are taken to be 2T-anti-periodic in the Euclidean time direction, 

X (x + 2T,x) = -x(x), x(x + 2T,x) = - X (x). (4.2) 

The orbifold reflection T about the Euclidean time xq = is defined by, 

T : x{x) -> (Tx)(x) = i7o75X(-^o,x), x{x) -> (xT)(x) = x(~x , x)?7o7 5 . 

(4-3) 

Following ref. |21| the gauge field is treated as an external field. It is first extended 
to the doubled time interval [— T, T], through 

U k (-x ,x) = U k (xo,x), U (-x - a,x) f = U (x), (4.4) 

and then 2T-periodically continued to all Euclidean times. Since T is an involution, 
one may define associated projectors and project on T-even and -odd fields as follows: 

X± = ~(l±T)x, X ± = x\(l±T). (4.5) 

Even and odd fields then satisfy homogeneous Dirichlet conditions at xq = 0, 

(1^707s)x±(0,x) =0, X±(0,x)(1^7 075 ) = 0. (4.6) 

Furthermore, due to the 2T-anti-periodicity in time, Dirichlet conditions with the 
complementary projectors are obtained at xq = T, 

(1 ± i7075)x±(T, x) = 0, x±(T, x)(l ± ^7075) = 0. (4.7) 
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T-even and odd anti-quark fields are defined analogously. For the orbifold projection 
to be consistent with the lattice theory, it is necessary that the reflection T commutes 
(or anti-commutes) with the Wilson-Dirac operator. The lattice action can then be 
decomposed according to 

Sf[x, X, U] = S f [x+ + X-,X+ + X-,U} = Sf[x+,X+, U] + Sf\x-,X-, U}. (4.8) 

To check whether this is indeed the case for the reflection (I4.3p . the crucial relations 
to note are 

TVq = -V T, TV = -VqT, (4.9) 

for the forward and backward covariant time derivatives. On the other hand, all spa- 
tial derivatives simply commute with T. Taking into account the 7-matrix structure 
in the Wilson-Dirac operator this then leads to [T, Dy/\ = 0, as required. 

As a consequence, one may consistently restrict attention to either the T-even or 
the T-odd fields. For instance, using only the T-odd fields, the fermionic functional 
integral is taken to be 

j D[x-,x-] eM-S f [x-,X-,U}). (4.10) 

Both T-even and T-odd fields are determined for negative times xo € [— T, 0], once 
specified for xo £ [0,T]. Hence, the field values at the negative times are no in- 
dependent degrees of freedom, and the integration measure over the odd fields in 
Eq. (I4.10p can be written as a measure over the fields at Euclidean times xq £ [0, T], 

D[x-,X-] = n d ^( x ) d ^( x ) II dX-(s)dx-(ar), (4.11) 

x 0<x <T 

where the non-Dirichlet components have been combined in the spinors </?(x) and 
<p(x), viz. 

p(x) = Q_X-(0,x) + Q +X -(T,x), (4.12) 
£(x) = X-(0,x)Q_+x-(Tx)Q+, (4.13) 

with Q± = ^(1 ±27075). Also the action can be reduced to the interval [0, T]: taking 
into account Eq. (|4.4p . the action for the fields at negative times xq € [— T, 0] turns 
out to be the same as for the fields at positive times. Eventually, the only trace left 
by the orbifold reflection is the structure of the Wilson-Dirac operator near the time 
boundaries. 

The case of a flavour doublet with boundary conditions given in terms of the 
projectors Q±, is naturally obtained by interpreting the even and odd fields as the 
components of a flavour doublet. More precisely, one sets 

i> = ^( X -\ V^ = v / 2(x-, X+). (4.14) 
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With this convention the fields tp and tp satisfy the boundary conditions (|3.5p . The 
prefactor has been included in order to obtain the correct normalisation of the 
kinetic term, once the action is reduced to the interval [0, T]. In the case at hand 
the dynamical degrees of freedom are the fields tp(x) and ip{x) for < xq < T as well 
as the non-Dirichlet components at the boundaries, i.e. Q-ip(0, x), Q + ip(T,x) and 
?/>(0,x)Q-, tp(T, x)Q+. The Dirac operator is implicitly defined when the action is 
reduced to the time interval [0, T]: 

Sf[x,X,U] = S f [x+,X+,U} + S f [x-,X-,U] 



1 4 

2 a 



E 



ip{x)(D w + m Q )tp(x) 



-T<x <T x 

° 4 E E tp{x)(V w + mo)ip(x). 

0<x <T x 



(4.15) 



More explicitly, writing Dyy as a time difference operator, 

aD w ifj(x) = -U (x)P^(x + aO) + Kip{x) - U (x - a6) t P+V'(x - a6), (4.16) 
with the time diagonal kernel being given by 



Ki/>(x) 



\ k=l 



tp(x) 



T C sw , CI 



3 

E 

fj,,u=0 



^F^ix^ix) 



(4.17) 



one finds 



\Q-KQ^{x) - Q_P^(x + aO) 
-P + Q^(x - aO) + Kip(x) - P^(x + aO) 
a£'vy' i / ; (2 ; ) 

-P+V(^ - aO) + JiTV(s) - P~Q+tp(x + aO) 
-Q+P+V(x - aO) + \Q + KQ + ^{x) 



if xo = 0, 

if xq = a, 

if a < xo < T — a, 

if xo = T — a, 

if x = T, 



(4.18) 

where the temporal links have been set to unit matrices for better readability. Note 
that this structure follows from the boundary conditions for both i/j and tp, and the 
explicit factor 1 /2 at the boundaries arises from the standard normalisation of the 
kinetic term. The orbifold construction also determines the Sheikholeslami-Wohlert 
term at the time boundaries. However, in practical applications it may be easier 
to set the coefficient of the Sheikholeslami-Wohlert term at the time boundaries to 
zero. This amounts to a change of a dimension 5 operator at the boundaries, i.e. a 
change in the action by a term of 0(a 2 ). 
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4.2 Orbifold reflection with an O(a) offset 

On a hypercubic lattice it is natural to distinguish between site and link reflections. 
The orbifold construction just described corresponds to using a site reflection about 
the time xq = 0. Orbifold constructions based on link reflections correspond to an 
offset by a/2 of the reflection point, i.e xq = ±a/2. Both options may be interesting 
in practice and are discussed in turn, beginning with link reflection about xq = —a/2. 
It is convenient to set T' = T + a. Then the starting point is the same as in 
Eq. (|4.ip . with the replacement T — >■ T", and the fermion fields are now taken to be 
anti-periodic with period 2T' . The orbifold reflection is defined by 

T+ : xO) -> «7o75X(-o - x , x), x(» -> x(~a - x , x)z7 7 5 . (4.19) 

and the external gauge field is extended to the interval [— T',T'], 

U k (-a - x ,x) = U k (x ,x), U (-2a - a; ,x) t = U (x), (4.20) 

and then 2T'-periodically continued. This implies that the spatial gauge fields at 
the time boundaries are duplicated and the time-like links connecting the doubled 
boundary layers are set to unit matrices. 

Again the orbifold reflection 7+ commutes with the Wilson-Dirac operator, and 
7+-even and -odd quark fields can be defined as previously, 

X± = ±(l±T+)x, x± = x|(l±T+). (4.21) 

However, it is obvious that these fields will not exactly satisfy the desired continuum 
boundary conditions. For instance, for the odd fields one finds, 

X-(0,x) = i7o75X(-a,x), X-( T , x ) = -«7o75X(^ + a, x), (4.22) 

i.e. boundary conditions would be obtained at xq = —a/2 or xo = T + a/2, where 
there are no lattice points available. One may however say that the boundary 
conditions are satisfied up to 0(a) effects. In fact, combining the even and odd 
fields in quark doublets ij) and t/> as in Eq. (|4.14p . one may write, in terms of the 
ordinary lattice forward and backward derivatives, 

Q+i 1 ~ ^0^)^(^)1*0=0 = 0. Q-i 1 + ^ a do)^{x)\ Xo =T = 0, 

$(x)Q + (l-±ad*)\ Xo=o = 0, iP(x)Q-(l + ^ad )\ Xo=T = 0. (4.23) 

Note that there is no problem with gauge invariance here, as the time-like links 
between the doubled boundary layer are set to unit matrices. 
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Once again the action and the dynamical variables can be reduced explicitly to 
the time interval [0, T], and the final result can be written as 



T 

Sf = a 4 ^2Y,^ x ^ V w + m )^(x), (4.24) 

xo=0 x 



where the Wilson-Dirac operator is specified by 



-U {x)P-ip{x + aO) + {K + i<y 5 T 6 P-)i)(x) if x = 0, 



aD\vip(x) if < xq < T, 

(K + ij 5 T 3 P + )i(j(x) - U (x - a6)^P + tp{x - 06) for x = T. 

(4.25) 

The dynamical field variables (i.e. the integration variables in the functional integral) 
are now given by the fields ip{x) and ip(x) for times < xq < T. Note also that the 
form of the Dirac operator T>\y can be obtained directly by acting with the infinite 
volume operator Dw on functions ip{x) for < xq < T, supplemented with the 
"syntactic extension" , 

ip(—a, x) = -z7o7 5 r 3 'i/'(0, x), ip(T + a, x) = ij -/ 5 T 3 ip(T, x), (4.26) 

and taken to vanish at all other times. Note that a similar statement can not be 
made for the orbifold construction based on the site-reflection T, as the expression 
of the Dirac operator in Eq. (14.180 also refers to the boundary conditions for ip. 

4.3 Orbifold construction with link reflection about xq = a/2 

We start again from Eq. (j4.1[) . but with the replacement T — > T' = T — a. The 
fermion fields are taken to be 2T'-anti-periodic and the orbifold reflection is defined 
by 

71 : x{x) -> «7o75X(a - x ,x), x{%) xifl ~ x)^7o75- (4.27) 

In this case, however, it is not possible to treat the SF gauge field consistently as an 
external field, as this would e.g. require to set 

U k (a,x) = £4(0, x), (4.28) 

i.e. the gauge field near the boundary would be equal to the gauge field at the 
boundary. However, this is not really important as the orbifold construction is no 
end in itself. It merely serves to teach us the structure of the Wilson-Dirac operator 
near the boundaries for the non-interacting theory. Hence the introduction of gauge 
fields can be postponed until the end, and whether or not the interacting theory 
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satisfies the boundary conditions will depend on the symmetries and the allowed 
counterterms. 

Proceeding in this way we obtain the action, 



T-a 




A ^2 ^2 i>(x)(D w + m )tp(x), 



(4.29) 



xo=a x 



where the Wilson-Dirac operator is specified by 



-U (x)P-tp(x + aO) + (K + ij 5 r 
aDw4>{x) 



= < 



if xq = a, 

if a < xq < T — a, 



(K + i-f 5 T 3 P + )^{x) - U (x - a6)tp + ^(x - aO) for x = T-a. 

(4.30) 

Note that the dynamical field variables in this case are ip{x) and i/j(x) at times 
< xo < T, which is analogous to the standard SF for Wilson quarks. 

4.4 The free quark propagator 

In order to check whether the boundary conditions are indeed satisfied in the free 
theory one may compute the free quark propagator in a time-momentum represen- 
tation, and compare to its expected continuum limit. There are various ways to 
proceed and some explicit results can be found in Appendix B. In fact the orbifold 
construction itself allows to calculate the propagator, taking the propagator on a 
hypertorus as starting point. The 3 different orbifold procedures can be combined 
by introducing the parameter r = ±1,0: we set T' = T + ra and identify the 
site reflection T with 7o- Performing the inverse Fourier transform in time of the 
propagator in four-momentum space, one obtains 




(4.31) 




(4.32) 



DwiP) = 



ip 1 + (m + \ap 2 ) 



(4.33) 



> 2 + (mo + \ap 2 ) 2 



Here, lattice momenta have been denoted by 



1 

p ll = - sin ap^ 




(4.34) 



17 



Note that there are IT' la allowed values for po, parameterised by an integer k, 

Po(k) = (2fc 2 + 1)7r , k = -T'/a, -T'/a + 1, . . . , T'/a, (4.35) 
so that the sum over po can be interpreted as a sum over k, 

T'/a-l 

E-E- E ( 4 - 36 ) 

PO po(fc) k=-T'/a 

All one needs to do to obtain the propagator in the orbifolded setting is to project on 
the appropriate 7^-even and -odd field components, and take the field normalisation, 
Eq. (|4.14p . into account. Combining the orbifold projection for both flavours by 
setting 

r T = diag(r T ,-r T ), (4.37) 



one obtains 



S;(x ,y ) = l(l-t T )Sl T '(x -ya)(l-%) (4.38) 
= Sl T '(xo-yo)-Sl T '(x + yo + Ta)i lol5 T 3 . (4.39) 

where the orbifold projection from the left (right) acts on 5p T (xq — yo) taken as a 
function of xq (yo)- Note that in going from Eq. (|4.38p to Eq. ()4.39p the behaviour 
under time reversal has been used, 

(t) = imSf H)*70T5- (4-40) 

It is now straightforward to check that the defining equations of the fermion propa- 
gator are satisfied, i.e. 

D^a 3 J>^ x -^(x ,y ) = a~% y , (4.41) 

X 

provided the arguments are restricted to < xo,yo < T for r = 0,+l and to 
< XQ,yo < T for r = —1. Furthermore, the homogeneous boundary conditions are 
indeed satisfied, for instance, 

Q+S;(0,y ) = Q + [sl T '(-y ) - Sf \-y - to)] , (4.42) 

vanishes for r = and is of 0(a) for r = dbl. Note that, formally, the boundary 
conditions are satisfied at xq = — ro/2, which may be used as a check, even though 
there is no lattice point at xq = ±a/2. 

A more explicit expression for the quark propagator can be obtained using the 
procedure of |40| . The resulting expressions are given in Appendix B and coincide 
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numerically with the above orbifolded expression^, provided one restricts the time 
arguments to < xq, yo < T' for r = 0, +1 and to < xq, yo < T' for r = —1. 

Finally, the orbifold procedure makes it clear that the free spectrum of j^t 1 !)^ 
is the same as obtained on a torus with 2T'-anti-periodic boundary conditions. More 
precisely, the eigenvalues are of the form 



A = ±^p 2 + (mo + lap 2 ) 2 . (4.43) 

In particular, the smallest (in modulus) eigenvalue is obtained for p = and is (for 
mo = and 8 = 0) given by 



IA • \- 2 ~ 
a 



( air 
sin — — 
V4T' 



~ 2T> ( 4 - 44 ) 



which converges to the known continuum result [2]. 
4.5 Orbifold construction, summary 

The orbifold construction yields the structure of the Wilson-Dirac operator near the 
time boundaries, such that the chirally rotated boundary conditions are realised in 
the free theory, possibly up to 0(a) effects. Using the link reflections 7±, the action 
of the Wilson-Dirac operator on the interval can be specified by a syntactic extension 
of the function space it acts upon. This is not the case for the site reflection T, where 
in addition one needs to refer to the boundary condition for the anti-quark fields. 
In all three cases the relation, 

( 75 T 1 P w ) t = 1$t 1 V w , (4.45) 

implies that the determinant of the Wilson Dirac operator is real. Moreover, a 
well defined eigenvalue problem is obtained for ^^Vy/- In the free theory, the 
eigenvalues are determined by the possible values of the time momentum component 
for anti-periodic boundary conditions on the doubled time interval. In particular, 
the lowest eigenvalue is determined by 7r/(2T"), where T' = T + ra and r = ±1,0 
for the three versions of the orbifold discussed above. 

Concerning mass terms, we have seen that the standard mass term can be added 
without any problem. Adding a twisted mass term i/i q 7sT 3 is, on the other hand, 
only possible by including a time profile following ref. [21]. Introducing the mass 
term in another flavour direction, e.g. i^ q 75T 2 is however possible without such a 
time profile, and corresponds to a twisted mass term in the standard SF basis, too. 



For r = 0, this is correct except for xo = yo = 0, T. However, this contact term in the 
propagator is never referred to provided the boundary quark and anti-quark fields are defined at 
xq = a, T — a (cf. Sect. 5). 
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In order to incorporate such a mass term in the orbifold construction it is best to 
start directly with flavour doublet fields and define the orbifold reflection directly 
for this case. The twisted mass term is then seen to commute with the spin-flavour 
matrix 27o7st 3 in the orbifold reflection and therefore with T T itself. 

5 Symmetries, renormalisation and O(a) improvement 

Given the fermionic lattice action with the correct boundary conditions at tree level, 
the next step consists in identifying the exact lattice symmetries and the allowed 
boundary counterterms. One then needs to choose a counterterm basis and discuss 
their lattice implementation and their effect on the basic fermionic 2-point functions. 

5.1 Counterterms of dimension 3 and 4 

The Ward identities, Eqs. (|3.31I3.32|) suggest that the symmetry breaking induced 
by the chirally rotated SF boundary conditions is quite similar to the breaking by a 
twisted mass term, except that the breaking terms are localised at the boundaries 
and come with dimensionless coefficients. Therefore it should not be too surprising 
that the symmetries of the lattice action in the chirally rotated Schrodinger func- 
tional are the same as in twisted mass lattice QCD [39], except for those space-time 
symmetries which mix time and spatial directions. In particular, there is charge 
conjugation, spatial lattice rotations, space and time reflections combined with a 
flavour exchange and global U(l) vector like rotations with generator r 3 /2. 

Given the symmetries, a list of fermionic operators of dimension 3 and 4 can 
be worked out. Their integrals over space, taken at xo = and xq = T define the 
possible counterterms to the lattice action that are needed to renormalise and 0(a) 
improve the boundary effects in the Schrodinger functional. We do not need to 
discuss the bulk counterterms here, as these are the same as in infinite volume |31j . 

Symmetrising with respect to charge conjugation we find the following operators 
allowed by the symmetries: there are 3 operators of mass dimension 3, 



01 

2 

03 



tpQ+tp, 



(5.1) 
(5.2) 
(5.3) 
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and further 8 operators of mass dimension 4: 

4 = 4>Q+"tkD k ^-i)D klk Q + ^, (5.4) 

_ ~ _ <— 

5 = ipQ-i k D k ip - ip D^kQ-ip, (5.5) 

6 = $Q +lQ D^-$D olQ Q + ^, (5.6) 

7 = ^Q-7oW-^B o7 oQ_V, (5.7) 

8 = ^Q+Do^ + ^SoQ+V?, (5.8) 

9 = ^Q_AW> + V'A)Q-V', (5.9) 

_ _ _ <— 

010 = ^Q+lQlkD k i> + ipD k -f k j Q + ip, (5.10) 

_ _ _ «— 

011 = ^Q-jo'jkDkip + tpD k -f k -f Q-^- (5.11) 

The dimension 4 operators are related by the equations of motion 

Jpi) = 0, $0 = 0, (5.12) 

which imply four relations 

04 + 06 = 0, (5.13) 

05 + 07 = 0, (5.14) 

8 + 0io = 0, (5.15) 

09 + 0ii = 0. (5.16) 

Furthermore, the combination, 

0io - 0ii = d k (i; 7 fc*75T 3 V0 , (5.17) 

yields a total spatial derivative and therefore does not contribute to the action. In 
summary, we have 5 relations among the 8 counterterms of mass dimension 4. No 

such simplification is possible for the dimension 3 operators. Hence one needs to 
consider 3 0(1) and 3 O(a) counterterms, which must be added to the SF lattice 
action with appropriately chosen coefficients. 



5.2 Formal continuum theory with inhomogeneous boundary conditions 

In order to understand how the boundary counterterms affect SF correlation func- 
tions one would ideally start from the lattice theory with inhomogeneous boundary 
conditions, following the steps taken for the standard SF |40j . Given the coun- 
terterms to the action, the renormalised and O(a) improved 2-point functions are 
obtained by differentiating with respect to the fermionic boundary values and bulk 
source fields, which are subsequently set to zero. Arbitrary SF correlation functions 
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can be built from these basic 2-point functions, and their renormalization and 0(a) 
improvement properties then follow, too. The main problem one faces is that the 
formulation of the lattice theory with inhomogeneous boundary conditions is not 
directly available. Indeed, our construction of the lattice theory based on the orb- 
ifold procedure yields homogeneous boundary conditions. While this is all that will 
be needed eventually, some alternative method is required to understand the role 
of the boundary counterterms. The approach taken here is based on the fact that 
the allowed counterterms of dimension 3 and 4 completely fix the structure of the 
unknown lattice theory to 0(a). This turns out to be sufficient to extract the desired 
information which may then be used to parametrize the renormalized lattice theory. 

As a first step let us go through the formal continuum theory with inhomogenous 
boundary conditions, specified by 

Q + ip(x)\ X0=0 = p(x), Q-tp(x)\x =T = p'(x), 

$(x)Q+\x =o = p(x), i)(x)Q-\ X0=T = p(x). (5.18) 

As in the case of the standard SF [2], the action is not just given by the bulk action 
but also requires local boundary terms. We thus make the ansatz 

Sf nt [^^,A} = J d 4 x^(x)i^(x) + J d 3 x [^(x)R ( X ^(x)] xo=0 

+ J d 3 x [ip(x)R T (x)iP(x)] xo=T , (5.19) 

with Rq(x) and i?r( x ) to be determined. Here, denotes the external gauge field, 
and the symmetrisation of the derivative implies that bulk and boundary parts of 
the action density must be separately invariant under charge conjugation. In order 
to determine the boundary terms we now follow the steps in ref. [2], and demand 
that the action must have smooth stationary points, i.e. classical solutions i/j c i and 
tp c \ which satisfy both the boundary conditions, Eqs. (|5.18|) . and the equations of 
motion, 

IP4> C \ = 0, $J/> = 0. (5.20) 
This is ensured if Rq and Rt are of the form 

flo(x) = -^ 75 r 3 + Q+r (x)Q + , R T (y) = - Z - 15 t 3 + Q_r T (x)Q_, (5.21) 

with arbitrary Tq and Tt, subject only to symmetry requirements. Imposing the 
continuum symmetries of the standard SF, in particular parity and flavour symme- 
tries, one finds that these terms must vanish, so that, 

i?o(x) = i2 r (x) = -^75T 3 . (5.22) 
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Not surprisingly, the resulting action coincides with the chirally rotated continuum 
action of the standard SF [2J. Using this action, we now consider the functional 
integral over the fermion fields, 

Z F [p,p;p',p';r ] ,fj} = [ Dty, $ exp {-Sf nt [^, A] + (fj, V) + , (5-23) 



which implicitly depends on the fermionic boundary values in Eq. (|5.18p . We have 
also included source fields 77 and fj in the bulk using the notation 



(77,^) = J d 4 xr](x)ip(x). 
The functional integration is best carried out after shifting variables, 

ip = i>A + v, i> = V>d + v, 



(5.24) 



(5.25) 



where the fluctuation fields v and v satisfy homogeneous boundary conditions, and 
the classical solutions of the field equations are explicitly given by 



d 3 y 



d 3 y 



s{x;0,yhoQ+p(y) - s(x;T,y)~f Q-p'(y) 

p'(.yhoQ+S(T,y;x) - p(y)7oQ_S'(0,y;x) 



(5.26) 
(5.27) 



Here, S(x, y) denotes the continuum quark propagator in the external gauge field 
A^. Integration over the v and v yields 



lnZ F = \nZ Ffi - S} wt [ip ch 4i, A] + (fj, S V ) + (fj, Vol) + (Ah v), (5-28) 



Zjr i0 = Z F [0,0;0,0;0,0], (Sr/)(x) = / d 4 y S(x, y) V (y). (5.29) 



with 



The dependence upon the source fields becomes completely explicit when insert- 
ing the classical solutions f|5.26|) and (|5.27p into equation (|5.28p . In particular the 
continuum action of the classical fields takes the form 



Sf^ ch ip ch A] 



d 3 xd 3 y 



p(*hoQ-S(0, x; 0, y)Q-7o/?(y) 

- p(x)7oQ-5 , (0,x;T,y)Q + 7op'(y) 

- p'(^hoQ+S(T,x;0,y)Q^ o p(y) 

+ p / (x)7oQ + 5(r,x;T,y)Q + 7op , (y) . (5.30) 
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It is convenient to follow ref. |40| and define fermionic expectation values through^ 

[0] F = \^OZ F \ . (5.31) 

Here the observables O may contain derivatives with respect to the source fields. 
Since the effective action is bilinear in the source fields, the fermionic expectation 
value of any observable can be given in terms of the basic fermionic 2-point functions, 
defined as all possible second derivatives of InZp with respect to the source fields. 
In order to be compatible with Eqs. (|3.19|) . for a = tt/2, one identifies 

5 5 
%(x)7oJ <H7o/»(x)J 

C'(x) o * C'(x) O 6 (5.32) 

<5[p'(x)7o] <H7oP'(x)] 



Using this notation, and, similarly, 



6 5 

il>{x) *+ prr, ip(x) O -jj-T-T, (5.33) 
077(3;) 077(0;) 

all the basic fermionic 2-point functions are given by 

(5.34) 
(5.35) 
(5.36) 
(5.37) 
(5.38) 
(5.39) 
(5.40) 
(5.41) 
(5.42) 

Thus, in the formal continuum theory, the insertion into Eqs. (|5.35I[574*2l) of the 
non-Dirichlet field components at the time boundaries (|3.19p . is indeed equivalent 
to their inclusion as functional derivatives, Eqs. (|5.32I5.33|) . 



[ip(x)Tp(y)} F 


= S(x,y), 




= S(x;0,y)Q_, 




= S(x;T,y)Q + , 


iaxMy))F 


= Q_5(0,x;y), 


[CW(y)]f 


= Q+S{T,x;y), 


[C(x)C(y)]F 


= Q_S(0,x;0,y)Q_. 


[C(x)C'(y)]F 


= Q_5(0,x;T,y)Q + 


[C'(x)C(y)]F 


= Q+S(T, x;0,y)Q_ 


[C'(x)C'(y)]F 


= Q + S(T,x;T,y)Q + 



5 The definition is such that the usual Euclidean expectation value is obtained by a subsequent 
integration over the gauge field, i.e. (O) = {[0]f)g, where the effective gauge field measure in (-)c 
includes both the pure gauge action and the fermionic determinant. 
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5.3 Including the counterterm action 

The correctly renormalised and 0(a) improved lattice theory is expected to converge 
to the formal continuum theory up to terms of 0(a 2 ). By adding the counterterms 
with some as yet undetermined coefficients one essentially defines Symanzik's ef- 
fective action for a generic lattice theory the formulation of which is not directly 
available. Obviously, the inclusion of dimension 3 boundary counterterms with ar- 
bitrary coefficients may alter the very continuum theory including the boundary 
conditions. This point will be addressed in Sect. 6, but we will ignore it here, as 
we are only looking for a structural result: by keeping arbitrary coefficients we may 
monitor their effect on the Dirac operator and the basic fermionic 2-point functions. 
This will then enable us to parameterize the counterterms in the lattice formulation. 

For definiteness, we use the equations of motion to obtain a counterterm basis 
for the massless theory, where all operators with temporal derivatives are eliminated. 
The resulting set of operators, 

01-5,010, ( 5 - 43 ) 
will be used at both time boundaries. The boundary counterterm action is given by 
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taking into account the behaviour under time reflections. Including the counterterm 
action, the total fermionic action becomes 



S 



total 
/ 



ocont 



+ Sffc + aSi. 



(5.51) 



The last term contains the volume O(a) counterterms [cf. Eq. 12.50 ]. Note that 
the elimination of counterterms containing time derivatives implies that O(a) effects 
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such as in Eqs. (|4.23p are absent. Hence, disregarding 0(a 2 ) effects, the boundary 
conditions can still be imposed at xq = and xq = T. In a first attempt we will 
thus try to maintain the boundary conditions (j5.18j) to 0(a). The first question to 
answer is whether there are still smooth classical solutions with the 0(a) enhanced 
action. Varying the action one immediately notices that there are three categories 
of counterterms. First, there are those which only refer to field components at 
the boundaries which are subject to the Dirichlet conditions, such as Q+ip and 
4>Q+ at xo = 0. This is the case for the counterterms oc 02,04. Then there are 
those which mix Dirichlet and non-Dirichlet components (oc c\,c\q) and finally the 
ones only referring to non-Dirichlet components (oc 03,05). Since the non-Dirichlet 
components are integration variables in the functional integral, the last category of 
counterterms can be absorbed in a re-definition of the Dirac operator. The same 
holds true for the bulk 0(a) improvement terms in S\. This leads to a redefined 
Dirac operator which we shall denote by Ipv> . Defining classical solutions, 

#< r ty d = 0, = 0, (5.52) 

subject to the boundary conditions (|5.18p . the question is whether the presence of 
the counterterms oc 01,02,04 and cio, allows for the classical solutions to be smooth. 
We need to check what kind of structures Rq and Rt correspond to the counterterms 
oc c\ , 02 , 04 and c\q , and whether these take the form given in Eq. (|5.2ip . This is 
indeed the case for the counterterms oc 02,04, which induce non- vanishing Tq and 
IV terms. However, the counterterms oc ci,cio cannot be cast in this form. Some 
thought reveals that one may take the counterterms oc c\ , c\q into consideration by 
a modification of the boundary conditions rather than by including them in the 
action. We set 

Q+V>(x)U =0 = Pr(x), Q_%l,( X )\ X0=T = p' r (x), 

$(x)Q+\x =o = Pr(x), tf;(x)Q-\ X0= T = p^(x), (5.53) 

where we have re-defined the boundary values as follows, 

p r (x) = v 7 ! - 2ci (1 - ci a7 fc Z) fe ) p(x), (5.54) 
p r (x) = p(x) (1 + c 10 aj k Dk) \/l - 2ci, cio = c 10 /(l - 2ci), (5.55) 

and analogously for the primed fields. When considered as a function of p, p, p' and 
p' the reduced action 

3 5 

s (r) = ^cont + J2 aSfl + c i aS fl + aS l- ( 5 - 56 ) 

i=2 i=4 
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coincides up to terms of 0(a 2 ) with the total action, Eq. (j5.51 j) . taken with the 
original boundary conditions. We re-define the Dirac operator as before, so that the 
fermionic action takes the form 

d 4 xi,(x)^ij(x)-± [ d 3 X [^(x)i 7 5T 3 ^(x)] 



\ / d 3 x [tp(x)i^ 5 T 3 ip{x)) 



x =T 



+ c 2 i>, A] + c 4 aSfilty, A}. 

The classical solutions ip c \,r an d tp c i r are defined by 

#< r tyd,r = 0, VH r ^ (r) = 0, 



(5.57) 



(5.58) 



subject to the modified boundary conditions Eq. (|5.53j) . Explicit expressions are 
similar to Eqs. (|5.26I5.2T|) . where the differences are accounted for by replacing the 
boundary values p—>p r etc., and the propagator S(x,y) — > S r (x,y), defined as the 
Green's function for Ip^ r \ Separating classical and fluctuation fields the fermionic 
action splits in two parts, 



S { p [Al,r + V, V? c l,r + V,A]= S ( j> [V, V, A] + Sj } [V>cl,n Al,r, A] 



(Or 



(5.59) 



The second term on the r.h.s. contains the whole dependence on the boundary source 
fields, and can be written in the form, 



c 2 j d 3 X L 
+ 2c 4 a J d 3 x 
+ / d 3 xd 3 y 



Pr(x)Q+Pr(x) + p',(x)Q_^(x) 

p r (x)Q + ^ k D k p r (x) + p' r (x)Q^~f k D k p' r (x.) 

p r (x)7oQ-S' r (0, x; 0, y)Q-7ofr(y) 

- Pr(x)7oQ-S' r (0, x; T, y)<3+7op'.(y) 

- Pr( x )7oQ+5'r(7 1 , x; 0, y)Q-70Pr(y) 

+ Pr( x )7oQ+5 , r(r,x;T,y)Q + 7op^(y) . (5.60) 



Introducing again the bulk source fields, we arrive at the expression of the fermionic 
generating functional, 



(r) 



InZ 



(r) 



' - SP [lpcl,r, 1pcl,r, A] + (ry, S r r]) 



+ (^,^cl,r) + (V'cl.r,??), 



(5.61) 



where the dependence on the source fields is completely explicit. Before taking 
derivatives with respect to the source fields one may partially integrate the terms 
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proportional to cio so that all the spatial derivatives act on the fermion propagator. 
It is then straightforward to compute all fermionic 2-point functions, e.g. 

[C(x)VJ(y)] F = v / r Tr 2^(l-cioa 7 fc£fc)Q-S r (0,x;y), ( 5 .62) 
[V>(x)C(y)b = VT^2^S r (x;Q,y)Q-(l + c w a^ k D k ). (5.63) 

By comparing with Eqs. (15.34I[5"74"2"1) . we thus obtain a clear picture of how the 
counterterms affect the basic 2-point functions and the Dirac operator. Based on 
this information one may now parametrize both the 0(a) improved Wilson-Dirac 
operator and the basic 2-point functions on the lattice. 



5.4 Lattice parameterisation 

We here focus primarily on the version of the xSF obtained with the orbifold reflec- 
tion T T with r = +1 (cf. Sect. 4). Hence, the fermionic fields at Euclidean times 
xq = and xq = T are integration variables in the functional integral and the coun- 
terterms in the chosen basis can be defined at xo = and xq = T. Starting with the 
Wilson-Dirac operator, the necessary counterterms are included by replacing T>w of 
Eq. (]4.25[) with V\y + SVw, where we define 

5V w 4>(x) = (5 x<h0 + S Xo>T ) [(«/-l) + (d.-l)oD,]^(x), (5.64) 

D s = I(V fc +V£) 7 fc- (5-65) 

On the lattice the fields are always subject to (a lattice version of) homogeneous 
boundary conditions (cf. subsection 4.2), so that ip and V> correspond to the fluctu- 
ation fields v and v in the continuum theory. This implies relations such as 

ijj(x)il>(x) =4>(x)Q-ip(x)\ +0(a) (5.66) 

a;o=0 lxo=0 

which have been used to simplify the parametrizationl (|5,64p . Note that 0(a) correc- 
tions to such continuum relations either re-define the coefficients of the existing 0(a) 
counterterms or modify 0(a 2 ) terms in the action, which are considered irrelevant 
in the present context. 

There is also considerable freedom in parametrizing the basic 2-point functions. 
We will here follow the convention used in the standard SF and obtain fermionic 
boundary fields from the fermion fields at times xq = a and xq = T — a, parallel 



6 When denning 8T>w for r = one needs to keep the projectors Q±, in order to project out 
the Dirichlet component of the fields at xo = and xo — T which are not integrated over in the 
functional integral (cf. Sect. 4). In the case r = —1, one may use the same parameterisation as in 
Eq. (|5.64p except that the Kronecker S's should be localised at xq = a and xo = T — a. 
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transported back to the boundaries at xo = and xq = T. In view of this, we 
introduce a shorthand notation for the parallel transported boundary-to-boundary 



propagators, 

S(0,x;0,y) = l7 (0 J x)5(a ) x;a,y)^o(0,y) t ) (5.67) 

S(0,x;T,y) = C/ Q (0,x)5(a,x;T - a,y)U (T - a,y) (5.68) 

5(r,x;0,y) = U (T - a, x?S(T - a, x; a, y)U (0, y)\ (5.69) 

S(T,x;T,y) = U (T - a, x)tS(T - a, x; T - a, y)U (T - a, y). (5.70) 

where the quark propagator is defined by 

(V w + 5V W + mo) S(x, y) = a~ A 5 Xty , (5-71) 

and we set niQ = m cr . The unrenormalised, 0(a) improved fermionic 2-point func- 
tions in the massless lattice theory are then parametrized as follows: 

Mx)$(y)] F = S(x,y), (5.72) 

bP(x)C(y))F = S(x;a,y)U (0, y y (l-d s ab s )Q-, (5.73) 

bP(x)('(y)}F = S(x;T-a,y)U Q (T-a,y)(l-d s ab s )Q+, (5.74) 

[C(x)$(v)]f = Q-(l + 4aD s )^ (0,x)5(a,x;y), (5.75) 

[C'(x)^(y)] F = Q + (l + d s aDs)U (T-a,^S(T-a,^,y), (5.76) 

[C(x)C(y)]F = Q-(l + d a aD.) [s(0,x;0,y) 

-(z f + d s aD s )<5(x - y)] (l - J s aD s ) Q_ , (5.77) 

[C(x)C'(y)]F = Q-(l + a>D s )S(0,x;T,y)(l-4aD s )Q + , (5.78) 

[C'(x)C(y)]F = Q+(l + 4aD s )5(T,x;0,y)(l-4aD s )Q_, (5.79) 
[C'(x)C'(y)]F = Q+ (1 + d.dD.) [S(T, x; T, y) 

-(z f + (f s aD s )<5(x - y)l (l - J s aD s ) Q+ , (5.80) 



where D s = ^ (V& + V£)7fc is the counterpart to D s in Eq. (|5.65|) . This defines 
the building blocks for any fermionic correlation function and therefore completes 
the set-up for the massless theory. It also implicitly defines the fermionic boundary 
fields as integrands in the functional integral. For instance, one has 

C(x) = (1 + J s aD s )C/ (0,x)Q_V(a,x), (5.81) 
C(x) = ^a,x)Q-U (x)^ (l-d s ah s ) , (5.82) 

except that these equations do not take into account the contact terms proportional 
to Zf and d s , which need to be subtracted whenever the possible Wick contractions 
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yield the 2-point functions [C( x )C(y)].F or [C'( x )C'(y)]F at coinciding spatial points 
x = y. 

The 2 renormalisation constants zt and Zf and the 3 0(a) improvement co- 
efficients d s , d s and d s implement lattice variants of the 5 counterterms Oi with 
i = 2, 3,4, 5, 10. The only missing counterterm is Oi, which can can be taken care 
of by a multiplicative renormalization of the fermionic boundary fields, 

(r = Z ( (, (r = Z c (, Cr = Z c C', Cn = Z ( (', (5.83) 

in complete analogy to the standard SF [3 31J. While the 0(a) improvement coeffi- 
cients d s and d s have their standard SF counterparts in ct and c s , the renormaliza- 
tion constants Zf and Zf and the coefficient d s are new and arise as a consequence of 
breaking the parity and flavour symmetries. However, the situation is not as bad as 
it seems, since correlation functions involving Zf and d s can be avoided in practice. 
Moreover, the counterterm proportional to d s is 7sT 1 -odd and only contributes 0(a 2 ) 
effects to 75T 1 -even correlation functions, thus rendering it irrelevant for most appli- 
cations. In practice, this leaves us with the tuning of the renormalisation constant 
Zf as the main new feature as compared to the standard SF. 

5.5 Away from the chiral limit 

The possible counterterm structures in the presence of mass parameters can be ob- 
tained by applying the equations of motion for massive fermions to the counterterm 
basis obtained for the massless theory. Possible counterterms are then given by those 
combinations which are invariant under charge conjugation. Considering the mass 
terms 

mtpip + Hqtpi'jsT 3 -^, (5.84) 
one finds that the boundary counterterms of dimension 4 have the form, 

{m, f i q }x{0 1 ,0 2 ,0 3 }, (5.85) 

and are thus proportional to the existing counterterms of dimension 3. They can be 
absorbed in a mass dependent rescaling of the corresponding renormalisation con- 
stants Z^, Zf and Zf. It should be emphasised that a non-zero standard mass term 
entails the loss of automatic 0(a) improvement. In order to maintain 0(a) improve- 
ment away from the chiral limit one thus needs to introduce all the usual 0(a) bulk 
counterterms, just like in the standard SF. This is to be contrasted with a twisted 
mass parameter, which preserves the 75T 1 -symmetry and is therefore compatible 
with automatic 0(a) improvement. 
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An interesting alternative mass term is obtained by introducing a twisted mass 
term in a different flavour direction^, 



n'^i^T 2 ^. (5.86) 

Such a mass term is invariant under the chiral rotation (|3.ip relating the standard 
and the chirally rotated SF, and is thus equivalent to a twisted mass parameter in 
the standard SF. Interestingly, in this case one needs a single countertem, which can 
be chosen to be 

(i> Q+HoIst 2 ^ ~ 4> Q-Hol5T 2 i?J . (5.87) 

Note that this term is 75T 1 -odd, so that it can be ignored, provided that m = 
and attention is restricted to 7sT 1 -even observables. Thus, no additional 0(a) 
mass counterterm is needed whilst automatic 0(a) improvement is maintained. The 
obvious drawback consists in the fact that the action ceases to be flavour diagonal. 

6 Renormalization and O(a) improvement 

The basic 2-point functions in the free theory can be used to determine the renor- 
malization constants 27 and zt and the 0(a) improvement coefficients d s , d s and 
d s at the tree level of perturbation theory. Some insight is provided for the role 
played by the renormalisation constant Zf. Automatic 0(a) improvement is then 
demonstrated with the free quark propagator in an Abelian background field. 

6.1 Interpretation of Zf and Zf 

The renormalization constants Zf and Zf are related to the counterterms O2 and O3 
which are 7sr 1 -odd. When chirally rotated back to the standard SF basis these 2 
operators become proportional to tjj' P±i'f5T 3 ^ 1 , which indicates a breaking of par- 
ity and flavour symmetry. Indeed, in terms of the standard SF fields the 75T 1 - 
transformation becomes a discrete SU(2) flavour transformation, 

ip' -> -tV $ -fir 2 . (6.1) 

Since parity and flavour symmetries are good continuum symmetries one expects 
that both renormalization constants are finite, i.e. scale independent functions of 
the bare coupling, with an expansion of the form 

z(g Q )=z^+g 2 Q z^ + 0(gt). (6-2) 
7 The author thanks Jenifer Gonzalez Lopez for reminding him of this option. 
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It is straightforward to determine their lowest order values, by comparing the basic 
2-point functions with a vanishing gauge background field to their expected contin- 
uum limit (cf. Appendix B). The free quark propagator on the lattice is obtained 
analytically for zt = 1. The fact that this yields the expected continuum quark 
propagator up to cutoff effects immediately implies 

zf = 1, (6.3) 

as the correct tree-level value. Furthermore, the absence of contact terms in the con- 
tinuum counterparts ()5.39|) and (|5.42p of the boundary-to-boundary 2-point func- 
tions (15.7715. 80p implies 

zf = \. (6.4) 

Note that these 2 finite renormalization constants restore the very 75T 1 -symmetry 
which enables automatic 0(a) improvement. A possible strategy to determine Zf 
and Zf therefore consists in the requirement that 2 suitable 7sr 1 -odd observables 
vanish exactly. Note that this strategy has implicitly been applied to obtain the 
results in Eqs. (|6.3|6.4p since the 7sr 1 -odd parts of the continuum 2-point functions 
all vanish identically. It should also be emphasized that Zf affects any observable 
via the Wilson-Dirac operator, whereas Zf only appears in the additive counterterm 
to the 2-point functions, Eqs. (|5.77|5.80p . Hence the determination of zj can be 
avoided if the observables of interest do not involve these 2-point functions. This 
can often be arranged by choosing appropriate quark flavour assignments [44j . 

Another perspective on Zf is obtained by considering an infinitesimal pertur- 
bation 5 around the lowest order value, Eq. (|6.3p . This would e.g. be the situation 
in perturbation theory with as the infinitesimal parameter. The first order cor- 
rection is equivalent to a single counterterm insertion in the correlation function of 
interest. Considering the quark propagator in the time- momentum representation 
and neglecting 0(a) artefacts, one obtains the expansion 

S U x °>y°)\ Zf =i +6 = S^ /2 \x o ,yo)-5\s^ /2 \x o ,0)S^ /2 \0,y Q ) 



+ si n/2 \x ,T)s£ /2 \T,y ) +0(a,o 2 ), (6.5) 



where the value of the angle a in the continuum propagator has been indicated by a 
superscript. Interestingly, the first order correction can be expressed as a derivative 
of the continuum propagator with respect to the angle a, 



5^ /2) (* , 0)5^ /2) (0, yo) + s£ ,2) (x , T)s£ /2) (T, y ) = -2-^S { p a) (x , yo] 



= 7r/2 

(6.6) 
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This gives rise to a complementary point of view regarding the role of zf. starting 
from a free theory with boundary conditions involving the projectors P±(a = vr/2) = 
Q±, the interactions tend to rotate the angle a away from tt/2. Since there is no 
lattice symmetry which protects the angle a = tt/2, one needs to tune Zf to maintain 
the angle at tt/2 all the way to the continuum limit. This interpretation is further 
corroborated by a numerical experiment: one may tune Zf such that the free lattice 
propagator Sp(xo,yo) approaches the continuum propagator defined at arbitrary 
values of a, as long as both time arguments are in the range < xq, yo < T. 

In conclusion, the finite renormalization constants Zf and Zf restore the 75T 1 - 
symmetry, which is interpreted as a flavour symmetry in the standard SF basis. With 
the chosen set-up for the fermionic boundary fields the tuning of Zf is equivalent to 
a renormalisation of the angle a such that the boundary conditions remain the ones 
with projectors Q±. In practice Zf and Zf can be fixed by tuning suitable 7sr 1 -odd 
observables to zero. In analogy with the 0(a) uncertainty in the definition of the 
chiral limit (cf. Sect. 2), one expects different determinations of Zf or Zf to differ 
by terms of 0(a). However, as in the case of the mass term, any 0(a) differences 
are expected to cause only 0(a 2 ) effects in 7sr 1 -even observables. 

6.2 O(a) improvement coefficients at tree- level 

To obtain 0(a) improvement at the tree level one needs to adapt the counterterm 
coefficients such that the continuum limit of the basic 2-point functions is approached 
with 0(a 2 ) corrections. All coefficients can be expanded in powers of g\ in analogy 
to the finite renormalization constants, Eq. (16. 2D . The first step should be the 
determination of d s since this counterterm enters the Wilson-Dirac operator and 
thus affects the fermion propagator in the bulk. Note that the exact expressions for 
the fermion propagator in Appendix B have been calculated without the counterterm 
5Vy/ (|5.64p . and therefore correspond to the choice d s = 1. Calculating the 0(a) 
effects in the propagator for this case yields, in terms of the coefficient functions of 
Eqs. (]B.31)B.32j) . 

TOj 

G±(p;x ,y ) = G±(p; x ,y ) + 12 {cosho;(xo - y ) 

2 cosh u)l 

=p sinho;Tsinha;(xo + yo — T)\ + 0(a 2 ), (6.7) 

H±(p;x ,yo) = H ± (p]x ,yo) + —^-—{smhuj(xo-yo) 

2 cosh loT 

=F sinhwTcoshwfco + Vo ~ T)} + 0(a 2 ). (6.8) 

We thus see that for r = the 0(a) effects vanish so that di°^ = 1 is the correct 
choice in this case. However, for r = ±1, there are uncancelled 0(a) effects, implying 
that = 1 is incorrect in this case. In order to retain d^ as a free parameter in 
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the lattice propagator to 0(a), one may proceed by calculating a single insertion of 
this counterterm into the propagator, viz. 

Sp(x ,y ) (0) = Sp(x ,yo) 

-a (4°) - l) [S { j /2 \x ,0) ip ^ k s£ /2) (0,y ) 



+ s£ /2 \x Q ,T)ip+>y k s£ /2) (T,y ) +0(a 5i ). (6.9) 



,2\ 



Hence the 0(a) effects in the first term of the r.h.s., explicitly given by Eqs. (16. 716. 8ft 
should be cancelled by the second term, which may be calculated directly in the 
continuum. After some algebra one finds that in all cases 

4°) =1-1 (6.10) 

appears to be the correct choic^l. This result has been confirmed by numerical 
inversion of the improved Wilson-Dirac operator in time-momentum space for a 
range of lattice resolutions, L/a. 

Having determined the 0(a) improved bulk propagator, it is straightforward to 
insert it into the basic 2-point functions. Again we limit ourselves to r = and 
r = +1. It turns out that a^ = in both cases. Finally, requiring the contact 
terms in Eqs. (|5.77I5.80|) to be absent at 0(a) yields 

4°) = 1. (6.11) 
6.3 A tree- level test of automatic O(a) improvement 

In the free theory the simplest non-trivial correlation function is the fermion prop- 
agator in a fixed external gauge field. We are thus led to consider "observables" of 
the type 

T r A = -a 3 J2 e ' ip{X ~ y) ^(y) r AQ±^(x)} F , (6.12) 

X 

where Ta denotes any of the 16 linearly independent Hermitian 4 x 4-matrices, 
7^, 1,75, i^^jBi Cfiv (for fj,,u = 0,1,2,3 and [x < v). Note that such bilocal quark- 
bilinear fields have a definite 7sr 1 -parity, which is is even (odd) if anti-commutes 
(commutes) with 75. Integrating over the fermions leads to 

I± A = tr{r A Q ± 5; =1 (x ,yo)} , (6.13) 

where the trace is over flavour, spin and colour indices, as these are contracted in 
Eq. (|6.12p . In the continuum limit, a chiral rotation relates these observables to 



8 Note that for this equation to hold with r = —1 it is assumed that the counterterm 5T>w (|5.64|l 
has been added to the Wilson-Dirac operator at Euclidean times xo — a, T — a. 
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standard SF correlators (cf. Sect. 3). With vanishing mass parameters and in the 
continuum limit one thus expects 

^r A ,cont =tr{T A P ± S { p a=0 \x ,y )}, if { 75 ,r A } = 0, (6.14) 

whereas 7sr 1 -odd observables vanish identically, due to flavour and parity symmetry 
(at a more technical level, in the standard SF basis, such observables are proportional 
to tr(T 3 ) = 0). It is instructive to consider the same observables in the standard SF 
at finite lattice spacing, which will be denoted by J' r . 

To take the continuum limit we keep all dimensionful parameters (xq, yo, T, . . .), 
fixed in units of L. Furthermore, we assume that boundary 0(a) improvement is 
correctly implemented, i.e. that the 0(a) improvement coefficients, d s or, in the 
case of the standard SF, c\, are set to their known tree- level values di°^ = 1 — r/2 
(|6.10p and c[°^ = 1 [31] . For T^'s which anti-commute with 75 one then expects 



1 r A 



„ . = 'ft.cont +«"»»<* (C»„ - 1)- + 0(0 2 ), (6.15) 

^.con t + 0(« 2 ). (6.16) 



{r\4,75}=o 



While the standard SF observables in Eq. (|6.15[) require c sw = 1 to be 0(a) improved, 
xSF observables in Eq. (|6.16p should reach the same continuum limit with a rate of 
0(a 2 ), regardless of the value of c sw - On the other hand, for 7sr 1 -odd observables 
one expects, 



1 r A 



1 A 



= 0, (6.17) 

[rA,7s]=o 

= const x (c sw - 1)- + 0(a 2 ). (6.18) 

[rA,7s]=o ■ L 



Like in the continuum, the exact flavour symmetry in the standard SF implies that 
75T 1 -odd observables are exactly zero [Eq. (|6.17p ]. While some of the xSF observ- 
ables vanish exactly, too, most are afflicted by 0(a) effects which can be reduced to 
0(a 2 ) by Symanzik improvement with c sw = 1, as indicated in Eq. (|6.18|) . 

In order to see the effect of the Sheikholeslami-Wohlert term at tree-level one 
needs the presence of a background gauge field. We will use the same constant 
Abelian SU(2) background gauge field which has been used for the definition of 
the SF coupling in the SU(2) Yang-Mills theory [41] . The fermion propagator in 
the continuum is known analytically in this case; however its lattice counterparts 
have so far only been calculated in the standard SF (cf. ref. [40]). We therefore 
resort to a numerical approach. We choose the set-up with r = +1, set T = L 
and p = 0, but keep a non-zero value 9 = 1.0. Specifying the Euclidean times 
xq = T/2 and yo = T/4 we choose the +- components with Ta = 70 and = 1 as 
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examples of 75T 1 -even and -odd quantities, respectively. Using the analytic result for 
the continuum propagator of Appendix B, their continuum limit can be calculated 
analytically, with the result, 



_ ^ cosh MT/2) - ^(r)] cosh^(T/4) 
^• cont " ^ cosh^(T) ' ' 



whereas -/"i~ cont = 0. The sum over colour components corresponds to the colour 
trace in the special case of the SU(2) Abelian background field. Computing these 
observables for c sw = and c sw = 1, and lattice sizes from L/a = 8 up to L/a = 
72, universality is indeed confirmed. In Figure 1 the relative deviations from the 
common continuum limit, rescaled by L/a, are plotted versus a/L. This extrapolates 
to a finite value in the case of 0(a) corrections, whereas a vanishing continuum 
limit indicates 0(a 2 ) or higher lattice artefacts. Finally, Fig. 2 shows the 7sr 1 -odd 
observable I±, multiplied by L/a. Hence, a finite value in the continuum limit 
indicates 0(a) effects. This is the generic situation unless c sw = 1 for which the 
observable becomes a pure lattice artefact of 0(a 2 ). In the free theory we have thus 
explicitly verified that automatic 0(a) improvement works as expected. 



7 Concluding remarks and outlook 

We have shown how the mechanism of automatic 0(a) improvement can be recon- 
ciled with SF boundary conditions. This means that any bulk 0(a) effects vanish in 
correlation functions which are even under the 75T 1 -symmetry. Any remaining 0(a) 
effects must arise from boundary effects and their elimination requires the tuning 
of a couple of boundary counterterms (proportional to d s and q), similar to the 
standard SF. In order to achieve this one needs to tune the coefficient Zf of a dimen- 
sion 3 boundary counterterm such that the physical parity symmetry is restored. 
One may say that the tuning problem for Zf is the price to be paid for avoiding 
the non-perturbative determination of c sw and of any 0(a) improvement coefficients 
for composite fields in the correlation functions (such as ca or cy for the flavour 
currents [3j). The chirally rotated SF provides a well-defined framework which can 
be seen as an alternative regularisation of the standard SF, with the advantage of 
reduced cutoff effects. This will benefit the calculation of step-scaling functions for 
which the continuum extrapolation will be made easier, and it will reduce the po- 
tential problem of introducing 0(a) effects into hadronic matrix elements through 
SF scheme renormalization factors. 

Since the formulation of the chirally rotated SF was given for a flavour doublet, 
a natural question to ask is whether it can be generalised to arbitrary flavour number 
N f . 
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Figure 1: The relative deviation from the continuum limit, rescaled by L/a, is plotted 
vs. a/L for the 7sT 1 -even observable (cf. text for further details). 



7.1 General N { 

Most of the discussion in this paper concerned the case of a quark flavour doublet, 
where a non-singlet axial rotation relates the correlation functions with those defined 
in the standard SF. In particular, invoking universality, one expects that the renor- 
malisation works out in the same way as in the standard SF, up to scale- independent 
(finite) counterterms, which are the consequence of breaking the continuum flavour 
and parity symmetry by the regularisation. This is further supported by the formu- 
lation with Ginsparg- Wilson quarks, where the same symmetries are restored a la 
Ginsparg- Wilson [25J. 

While the inclusion of further doublets and thus the formulation for even flavour 
numbers is straightforward, the situation is a bit more complicated for odd Nf. First 
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Figure 2: The 0(a) effect in the 75T 1 -odd observable if with c sw = is rendered 
visible by rescaling the observable with L/a. The effect of 0(a) improvement with 
c<= w = 1 reduces the lattice artefacts to 0(a 2 ). 



of all, the discussion of the symmetries, and hence the counterterm structure, relied 
on a flavour exchange between up- and down-type flavours combined with space- 
time reflections. To apply this analysis to the single flavour case one may use a 
spurion field in the Wilson-Dirac operator, which switches sign under a space or 
time reflection, and thus takes over the role of the flavour exchange symmetry. The 
counterterm structure found in this way would then be the same as in the doublet 
case, appropriately reduced to either the up- or down-type flavour. This may in fact 
be generalised to any number Nf of flavours, with all flavours being of the up-flavour 
type, for example. As the relation to the standard SF is lost in this case, one would 
then tend to describe the symmetries as an exact SU(iVf ) flavour symmetry, with all 
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axial symmetries being broken by a singlet twisted mass term with 0(1) coefficient, 
localised at the boundaries. In addition to the axial symmetries also parity and 
time reflection symmetries are broken explicitly by the boundary conditions. Note, 
however, that the interpretation that Zf restores a continuum symmetry is lost in 
this case, and it is unclear to the author whether Zf remains scale-independent in 
this situation. A further, more practical question is whether such a set-up would be 
suitable for numerical simulations. 



7.2 Numerical simulations including the fermion determinant 

Currently, most numerical simulations are based on some variant of the Hybrid 
Monte Carlo algorithm [42], and the fermionic determinant is thus required to be 
real and positive. For the Wilson-Dirac operator in the two-flavour case with xSF 
boundary conditions this property can easily be established. The hermiticity of 
75T 1 Pv^; implies for its flavour components, 



2V = diag Xy,X>. 



,(1) t,(2) 



'W 



75^75 



V 



(7.1) 



and therefore, 

det (Vy/ + mo) 



det ( V$ + m ) det (v { $ + m 



det 



pW+mo) [V#+mo 



> 0. 



(7.2) 



which has exactly the structure required for the standard algorithms. Twisted mass 
terms can be added, too. Allowing for two kinds of twisted mass terms, 



^q75T 3 +^q75T 2 , 



(7.3) 

one realizes that the Wilson-Dirac operator is still 75T 1 -hermitian so that its de- 
terminant must be real. Reducing the determinant in flavour space the r.h.s. of 
Eq. (|7.2j) generalises as follows, 



det 



1>w + m o 



(2) 



(7.4) 



Note that /i q multiplies a hermitian operator localised at the boundaries, which is 
diagonal in space-time and gauge field independent, at least for the set-up with 
r = ±1 (cf. Sect. 4). 

The situation is less clear with odd numbers of flavours. In the free theory 
the determinant for a single flavour turns out to be real, even in the presence of a 
spatially constant Abelian background field. However, a numerical check on a small 
lattice with a random gauge configuration reveals that the determinant is in general 
comple^. Hence, some additional idea is needed in order to enable simulations of 



I thank B. Leder for performing this numerical check using MATLAB. 
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odd flavour numbers with ^SF-type boundary conditions. A simple alternative in 
this case would be a mixed formulation, where an even number of quarks is repre- 
sented by chirally rotated quark doublets, whereas a single quark is implemented 
with standard SF boundary conditions. In the continuum limit this formulation will 
be equivalent to the standard SF with an odd number of flavours. Moreover, auto- 
matic 0(a) improvement would partially apply: while the Sheikholeslami-Wohlert 
term in the action would be required, correlation functions involving only chirally ro- 
tated quark flavours would be automatically 0(a) improved up to boundary effects. 
The price to pay is an additional 0(a) improvement coefficient at the boundaries 
since now both d s and ct would be required. 

7.3 Outlook 

The chirally rotated SF offers a number of applications. Automatic 0(a) improve- 
ment will be particularly welcome when applying SF schemes e.g. to four-quark 
operators for which a non-perturbative determination of the 0(a) counterterms 
seems impractical. Furthermore, in the standard SF, the 0(a) improvement by 
the Sheikholeslami-Wohlert term often cancels large cutoff effects. In order to apply 
SF schemes to operators used in twisted mass QCD without the Sheikholeslami- 
Wohlert term (as currently used by the European twisted mass collaboration [43J), 
the chirally rotated SF is a natural choice, and may indeed be the only practical 
one. Further applications include multi-flavour QCD-like theories or lattice gauge 
theories with different gauge groups and fermion representations, where the deter- 
mination of the Sheikholeslami-Wohlert coefficient c sw can be avoided. Moreover, 
one may expect a reduced sensitivity to the precise definition of the critical quark 
mass. Finally, the chirally rotated SF provides new methods to determine finite 
renormalisation constants such as the axial current normalisation constant Z\, and 
some of the usual 0(a) improvement coefficients [44| . Since standard Symanzik 0(a) 
improvement of Wilson fermions remains relevant for many collaborations, any im- 
provement over the traditional methods based on Ward identities may indeed be 
welcome. 

However, before these applications can be envisaged, one would like to per- 
form a thorough non-perturbative test of the chirally rotated SF. In particular one 
needs to acquire some experience regarding the tuning of zt and m cr |45|44j . One 
would like to verify that the boundary conditions are correctly implemented in the 
interacting theory and that automatic 0(a) improvement works out as anticipated. 
Furthermore, universality could be checked by comparing with correlation functions 
defined in the standard SF. A non-perturbative study addressing all these questions 
is currently being performed in quenched lattice QCD [33]. Further checks will be 
carried out in perturbation theory [36]. In particular, one would like to obtain per- 
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turbative estimates for d s and ct, as these improvement coefficients are difficult to 
determine otherwise. 
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A Absence of zero modes 

The absence of zero modes for the massless continuum Dirac operator in a smooth 
external gauge field An(x), can be established as in the case of the SF with standard 
boundary conditions [24] . We first deal with the case of a single quark flavour with 
Q± boundary conditions, and then consider the case of a quark doublet at generic 
values of the angle a. 

A.l Single flavour with Q± boundary conditions 

Consider the massless Dirac operator If) = j^D^ + A^) for a single flavour, acting 
in the Hilbert space of square integrable functions <p(x) on a hyper cylinder, subject 
to the boundary conditions 

Q+^(0,x) =0 = Q_(^(T,x), Q± = ^(l±i7075)- (A.l) 
With respect to the natural scalar product, 

(cp, X ) = J dx j d 3 xv?(x)t x (x), (A.2) 

the operator ilp is hermitian. A well-defined eigenvalue problem is thus obtained, 

ilpw n = \ n w n , (A.3) 

with real eigenvalues A n and smooth eigenfunctions w n , which span the Hilbert 
space. In the following it is assumed that the eigenfunctions have been chosen such 
as to form an orthonormal basis of the Hilbert space, i.e. 

(w n ,w m ) = S nm , ^w n (x)w n (y) ] = <5 (4) (x - y). (A.4) 

n 
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Next, consider the classical solution of the Dirac equation with inhomogeneous 
boundary conditions, 

#>Vci(x) = 0, Q+Vcl(0,x) =p(x), Q_Vci(T,x) = p'(x). (A.5) 

In order to obtain an explicit expression for the classical solution, one needs the 
propagator in the orthogonal subspace of the zero modes, 

S" L (x,y) = -i £ "" ( ^ (y)t . (A.6) 

n:A n ^0 

Assume that there are a finite number of iVo zero modes, w ni (x), i = 1, . . . , Nq. The 
propagator in the orthogonal subspace then satisfies the equation 

N 

lpS ± (x,y) = 5^(x - y) - £ w ni {x)w ni {y)K (A.7) 

i=l 

Furthermore, since any finite superposition of zero modes must be smooth, the 
discontinuity for S ± (x,y) is given by 

S>(x ,x,x ,y) - S<(x Q ,x,x ,y) = 7 5 (3) (x - y), (A.8) 

where the subscripts refer to the ordering of the time arguments in the propagator. 
Using these properties, the classical solution can be written as follows, 

^ci(x) = J d 3 y [S x (x,0, yhoQ+p(y) ~ S ± (x,T,y) l0 Q_p'(y)\ , (A.9) 

showing that it is linearly related to the boundary values. Finally, note that any 
of the zero modes, w ni (x), is a classical solution of the Dirac equation, and can be 
written in the form (|A.9p . Since its boundary values vanish, the zero modes must 
vanish, too, i.e. w Ui = 0, i = 1, . . . , Nq. 

The case of opposite boundary projectors is completely analogous. Hence, the 
absence of zero modes is guaranteed with any number of quark flavours, and for any 
mixture of up and down type flavours. 

A. 2 Flavour doublet at generic values of a 

Essentially the same argument works for a flavour doublet at any value of a, except 
that the mathematical structure is a bit more complicated. The boundary conditions 
for ip and at fixed, but arbitrary a give rise to the definition of two different Hilbert 
spaces, 

U± = {(f> : P±(±a)<t>\ xo=0 = = P T (±a)cf>\ X0 =T}, (A.10) 
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and both Jp and 75T 1 can be seen as mappings between these two spaces. Hence, 
75T 1 ^Z) is a hermitian operator in either space, with a complete set of smooth eigen- 
functions and real eigenvalues. Denoting its orthonormal eigenfunctions in % + by 

7 5 T 1 ^V n = X n ip n , (A.ll) 

it follows that the functions (p n = ^^r l tp n form an orthonormal basis in and 

Tp"ihT X <Pn = IpVn = A„7 5 tV„ = \ n <p n . (A. 12) 

Completeness in either Hilbert space then implies 

S 4 ( x -y) = c Pn(x)ip n (y)^ = ^2 ^WifttW 1 '■ (A.13) 
n n 

To relate the classical solution of the Dirac equation, Tp ipd(x) = 0, to its boundary 
values, 

P+(a)^ d (0,x) = p(x), P_(a)^ d (T,x) = p'(x), (A.14) 
one needs again the inverse of If) on the non-zero modes, 

sH*,y)= E (A.15) 

Assuming iVo zero modes, this propagator satisfies 

N 

l/)S ± (x,y) = 5^(x - y) -J^MxtfM 1 . (A.16) 

i=l 

The discontinuity is again described by Eq. (|A.8j) . and the classical solution takes 
the form 

Al(x) = I d 3 y ^ ± (x,0,y) 70 P + (a)p(y) - S ± (x,T,y) l0 P^(a)p' (y)] . (A.17) 



Again, any of the zero modes solves the Dirac equation and satisfies homogeneous 
boundary conditions, and must therefore vanish identically. Note that this argument 
also holds at a = 0. In particular, the flavour permutation r 1 is not required in this 
case, so that the proof applies virtually unchanged to the case of a single quark 
flavour with standard SF boundary conditions. 



B The free quark propagator 

We collect a few formulae for the free propagator both in the continuum and on the 
lattice. For the sake of notational simplicity the propagator is generically denoted 
by S and any additional parameter dependence or indeed whether it is the lattice or 
continuum propagator will be either clear from the context or otherwise indicated. 
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B.l Continuum quark propagator for generic values of a 

The free quark propagator in the continuum theory can be easily obtained using 
standard methods. Allowing for standard and twisted mass terms, the quark prop- 
agator satisfies, 

(jp + m + ^ q75 r 3 ) S(x, y) = 5 (4) (x-y), (B.l) 
and the boundary conditions, 

P + (a)S(x,y)\ X0=0 = 0, P-(a)S(x,y)\ Xo=T = 0, 

S(x,y)P-(-a)\ yo=o = 0, S(x, y)P + (-a)\ yo=T = 0. (B.2) 

The boundary conditions in the second argument follow e.g. from the fact that 
the Dirac operator is 75T 1 -hermitian (cf. Appendix A). For the same reason the 
complementary components satisfy (modified) Neumann conditions, i.e. 

(D - M a )P_(a)S{x,y)\ X0=0 = 0, (D + M a )P + (a)S(x,y)\ Xo=T = 0, 
S(x,y)P+(-a){D o -M a )\ yo=0 = 0, S(x,y)P-(-a)(D + M a )\ yo=T = 0, (B.3) 

with 

M a = m cos a + /i q sin a. (B.4) 

These equations hold for arbitrary smooth gauge field A^x). Setting the gauge field 
to zero, we use L-periodic spatial boundary conditions and include the phase 9 as 
an Abelian background field in the spatial derivative: 

J/> = Oojo + (dk + i^j 7fc . (B.5) 

Spatial translation invariance suggests to pass to the time momentum representation, 

S(x, y) = L~ 3 e ip{x ' y) S p (x ,y ), (B.6) 
p 

where the (infinite) sum is over all allowed momenta of the form p = 2-rrn/L, with 
n being a triple of integers. In time- momentum representation one then has 

(do - -4) S p (x , y ) = 70^0 - yo), (B.7) 

with 

^=(^7^-^ + ^75^)70 = ^, pt=Pk + 0/L. (B.8) 
There are many ways to represent the free propagator. A very compact form is 

S P (xo,yo) = e XoA {—^-P-(a)e TA + e(x - yo )}e-y° A 10 , (B.9) 
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where 9(t) denotes the Heaviside step function, and 

Ma 

Using A 2 = (J 1 one may easily evaluate the exponentials and obtain a more explicit 
representation, 



n(T) = cosh(wT) + — sinh(uT), u = ^/(p+) 2 + m 2 + n\. (B.10) 



Sp(x ,yo) = j P-(-cx)H + + ^P + (-a)H- 

-(A - M a ) 7o P_(-a)G+ - (A + M a ) 7o P + (-a)G_, (B.ll) 

with coefficient functions G± and i/± defined by 

G+(p;x ,y ) = ^-i r {(M a - W )(e^-^l- T )-e^' + ^- T )) 

+ (M a + w) ^e^^ ^ !-^ - e-^o+w-r)) J, (B.12) 

G_(p;x ,yo) = G + (p; T - x , T - y ), (B.13) 

-^±(p;^o,yo) = (-<% T M a )G±(p;x ,y ), (B.14) 

and 

i? a = (M a + w)e wT - (M a - w)e _wT . (B.15) 

Note that the coefficient functions are indexed such as to conform with the projectors 
for a = 7r/2, where P±(— n/2) = Q T . They enjoy the following properties: 

G±{p;x ,y ) = G±(p;y ,x ), H±(p;x ,y ) = -H T (p;y ,x ). (B.16) 

The Dirichlet boundary conditions imply 

G + (p;0,y ) = = G_(p; T, y ), (B.17) 

and 

J ff_(p;0,y ) = = J ff + (p;yo,0), H + (p; T, y ) = = tf_(p; Vo, T). (B.18) 
Of particular interest is the massless case, for which the formulae simplify, 

G ± (p;x ,y ) = — ^_ {-ainh^xo - »b| -T) 

±sinhu;(xo + yo-r)}, (B.19) 

H±(p;x ,y ) = 2 cosh lot { £ ( X ° ~ ^o) cosha;(|xo - yo\ -T) 

Tcoshw(x + y -T)}, (B.20) 

where e(t) denotes the sign function. 
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B.2 The massless continuum quark propagator in an Abelian and spa- 
tially constant background field 

In this subsection, an explicit expression is given for the massless continuum quark 
propagator in an Abelian SU(iV) background field B^. The latter is taken to be 
spatially constant, i.e. = B^(xq), and in the Bq = gauge. Spatial translation 
invariance then allows to pass to the time-momentum representation, where the 
propagator satisfies, 

(<9o7o + [ipt + B k (x )] 7 fc ) Sp(x ,y ) = S(x - y ). (B.21) 

We are interested in the background fields used for the running coupling |l)41)47j or 
the determination of the Sheikholeslami-Wohlert coefficient |48| . These are of the 
form 

B k (x ) = C k + ^{C' k -C k ), (B.22) 
where the boundary gauge fields are taken to be 

C k = i<j>/L, C' k = i<j>'/L, (B.23) 

with diagonal N x iV-matrices <f> and eft' in colour space. Since these are chosen to 
be independent of the spatial index k, it is convenient to restrict attention to the 
special case where also the momentum components are k- independent, i.e. p~£ = p + 
for k = 1,2,3. Then, defining 

S = -j= (71 +72 + 73 ) 70, (B.24) 

«(t) = *^(p + + i + 2lr^' (R25) 

the result for the quark propagator can be written in the form 

S P (x ,yo) = e<^ o)s {-^^(T/'^^ + 9{x ° " Vo) } ( B - 26 ) 

where 9(t) denotes the Heaviside step function. Using S 2 = 1, one may easily obtain 
the equivalent expression 

S p (x ,y ) = 7 P-(-a)i?+ + 7 i 3 +(-a)-i^- 

/(T)S 7o + P + (-a)G_} . (B.27) 
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Here, to' if) = duj(t)/dt, and the coefficient functions are defined by 
Xo,yo) — 2a;' (T) cosh u(T)l ^ ~*~ U ^° ^ 



sinh [e(x - yo) {w(a;o) - ^(yo)} - w(T)] } , (B.28) 



#±(p;^o,yo) = ^^-^yjTcosh [w(x ) + w(y ) - ^(^)] 

+e(x - yo) cosh [e(i - y ) {w(x ) - w(y )} - w(T)]j(B.29) 

Here all functions are diagonal in colour space, with colour components being ob- 
tained by specifying the colour components fa and <\>\ and thus Wj, with i = 1, . . . , N. 
Setting fa = fa = amounts to the replacement of Ui(t) — > oj(p + )t and these ex- 
pressions then reduce to the ones given in Eqs. (1RT91) and (lB~20i) . for the special 
case with pf = pj = pg = p + . 

B.3 The free quark propagator on the lattice 

A straightforward way to calculate the free quark propagator on the lattice proceeds 
via the classical solution of the lattice Dirac equation, following [3D], where the 
inhomogeneous boundary conditions are imposed either at Euclidean times xq = 
0, T, or appropriately off set by ±a/2. After spatial Fourier transform, the result 
can be written in the form 

Sp(x , yo) = 2A{p+) Sxo,yo + 10 ® +H + + ^°Q- H - 

+ {-Wtlk + M(p+)) [Q + G r + + , (B.30) 

where the coefficient functions are given by 

Gi(p;z ,yo) = 2^+)^^ {~ *"MI*o - yol - ?) 

±sinhw(x +y -T)\, (B.31) 

#±(p;zo,yo) = 2 yi(p+)cosh^r / {^ a;o ~ ^o) cosnuj (l x o ~ yo\ ~ T> ) 

=Fcoshcj(x + yo-r)}. (B.32) 

Here, T' = T + ra denotes the effective time extent for the fermionic degrees of 
freedom, for the orbifold construction with r = and r = ±1 (cf. Sect. 4). With 

A(p) = l + amo + ^p 2 , (B.33) 
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the one-particle energy oj is defined by 



sinh 




{ 



q 2 p 2 + {A(p) - l} 2 \ * 
4A(p) J 



(B.34) 



The notation is further specified by 



aM (p) 



— sinh(aw), 
a 

A(p) — cosh[aa;(p)], 



(B.35) 
(B.36) 



and the above lattice version of the sign function, e(t), is given by 



1 



if t >0 



e(t) = { 







if t = 



(B.37) 



1 if f < 0. 



These expressions for the free quark propagator have been checked numerically by in- 
verting the corresponding Dirac operators at fixed momentum, and semi-analytically 
by comparing with Eq. f|4.39j) . For r = a difference with Eq. f|4.39j) is found for the 
contact terms at the boundaries, i.e. for xq = yo = and xo = yo = T. While the 
difference can be easily worked out, the result is not needed here. In fact, these con- 
tact terms are never referred to provided the fermionic boundary fields are defined 
as in Sect. 5. 
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